Anomalous high energy dispersion in photoemission spectra from insulating cuprates 
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Angle resolved photoelectron spectroscopic measurements have been performed on an insulating 
cuprate Ca2Cu02Cl2- High resolution data taken along the V to (7T,7r) cut show an additional 
dispersive feature that merges with the known dispersion of the lowest binding energy feature, 
which follows the usual strongly renormalized dispersion of « 0.35 eV. This higher energy part 
reveals a dispersion that is very close to the unrenormalized band predicted by band theory. A 
transfer of spectral weight from the low energy feature to the high energy feature is observed as the 
T point is approached. By comparing with theoretical calculations the high energy feature observed 
here demonstrates that the incoherent portion of the spectral function has significant structure in 
momentum space due to the presence of various energy scales. 

PACS numbers: 71.27.+a, 78.20.Bh, 79.60.Bm, 74.72.-h 



Since the early days of high temperature supercon- 
ductivity (HTSC), it has been suggested, based on the 
cuprate phase diagram, that the microscopic origin of 
HTSC may hinge on how charges move in a background 
of strong antiferromagnetic (AF) interactions. Thus, 
a natural starting place is to learn how a single hole 
behaves in a half-filled AF system. This is precisely 
what is measured by angle resolved photoelectron spec- 
troscopy (ARPES) on the half-filled cuprates. Indeed, 
such considerations have motivated many ARPES studies 
on parent HTSC compounds such as Sr2Cu0 2 Cl 2 j 1 i 2 i 3 i 4 
Ca 2 Cu0 2 Cl 2 ^ Nd 2 Cu0 4 £ and lightly doped La 2 Cu0 4 . 7 

While Local Density Approximation (LDA) calcula- 
tions predict these systems to be metallic with the 
width of the lowest filled energy band ~2 eV^ optical 
measurements indicate they are insulators with a band 
gap of ~1.5 eV(2iiS which results from the strong on- 
site Coulomb repulsion. ARPES results have revealed 
many other important features of insulating cuprate 
compounds. Energy Distribution Curve (EDC) analy- 
sis has revealed a broad low-energy charge transfer band 
feature, 1 which is commonly referred to as a Zhang-Rice 
singlet. It has a band width of 0.35 eV which has been 
related to a renormalization by the exchange interaction 
J. The dispersion has been well described by the t-t'-t"- 
J models In addition, results from Ca 2 Cu0 2 Cl 2 show 
a rf-wave-like dispersion and a feature in the occupation 
probability n(k) that was termed as a 'remnant' Fermi 
surfaced 

All these results were obtained with a lower angular 
(that is, momentum) resolution of 2° compared to what 
is commonly used now (0.25°), yet made important con- 
tributions in understanding hole motion in the AF back- 
ground. As has been found across the phase diagram of 
the cuprates, the dramatic increase in photoemission res- 
olution has revealed new features in the data which were 
previously unresolved^ 1 - In this paper, we present high 



k- resolution data from Ca 2 Cu0 2 Cl 2 which reveal yet an- 
other surprising result. The lowest energy excitations 
near (ir/2, 7r/2) remain consistent with the above descrip- 
tion, but at higher binding energies we observe a broad 
rapidly dispersing feature which tracks the unrenormal- 
ized band dispersion. We note that there is significant 
transfer of spectral weight from the low energy feature to 
the high energy one as the T point is approached. This 
description is consistent with several different numeri- 
cal calculations r^^^Sii 2 ^ and our data experimentally 
confirms this high energy electronic structure. We will 
discuss our observations in the context of the current 
theoretical understanding. 

Ca 2 Cu0 2 Cl 2 single crystals were grown by a flux 
method. 12 Experiments were performed at beamline V of 
the Stanford Synchrotron Radiation Laboratory (SSRL). 
Samples, which were oriented ex situ prior to the exper- 
iments by Laue x-ray diffraction were cleaved in situ at 
a base pressure better than 5 x 10~ n torr. Spectra were 
taken with 25.2, 22.4, and 15.5 eV photons, the total en- 
ergy and angular resolution was typically 30 meV and 
0.25°, respectively, and no sign of charging was observed 
in these samples. 

Figure 1(a) shows EDCs of the lowest energy feature 
along the nodal direction (that is, k x = k y ) from tetrag- 
onal Ca 2 Cu0 2 Cl 2 . It is seen that a broad peak develops 
near (7r/4, tt/4), gains spectral weight, and disperses to- 
ward the lower binding energy side as k increases from 
r to (7r/2,7r/2). It then suddenly loses intensity near 
(tt/2, 7r/2). Furthermore, it leaves no spectral intensity 
at the Fermi energy Ep, as expected for an insulator. 
The total dispersion of the peak is about 0.35 eV. Note 
also that there is very little spectral intensity near the 
r point. These are the salient features of the insulating 
parent HTSC compounds and have been discussed in the 
earlier literature^ 1 " 

Our high fc-resolution data, however, reveal a feature 



2 




2.0 1.5 1.0 0.5 0.0 0.2 0.4 0.6 0.8 1 



Binding Energy (eV) Momentum (0,0)-» (ji.ji) 



> 

CD 



CD 
C 

C 

in 




(71/2,71/2) (0,0) (7t/2,7t/2) 

Momentum (k x =k y ) 



FIG. 2: Intensity plot of data shown in figure I as functions 
of the binding energy and momentum. The data was sym- 
metrized around the F point. Also shown on the plot are 
the dispersions obtained by following the peak positions of 
the MDCs (solid line) and the EDCs (circles and triangles). 
The results are compared with the shifted dispersion from the 
LDA calculation (dashed line). 



FIG. 1: ARPES data from Ca 2 Cu0 2 Cl2 along the r-(7r,7r) 
direction, with 25.2 eV photons, (a) EDC's with the k point 
for each curve spaced evenly along the r-(7r,7r) line, (b) MDC 
curves of the ARPES data shown at equal binding energy 
intervals, from 2 eV (bottom) to the chemical potential (top). 
The large spectral weight near (n,n) for high binding energy 
MDCs is due to the 02 P7r states in the main valence band. 



that was not previously observed. The details can be 
seen better with an intensity plot. The data in figure 1 
are symmetrized about the T point and plotted in Figure 
2. The low energy feature described above is clearly seen 
between k — (n/4, 7r/4) and (tt/2,tt/2) with total dis- 
persion of about 0.35 eV. Examining the image, one can 
clearly see that there is a fast dispersing feature above 
0.8 eV binding energy, which is in sharp contrast with 
the known and slowly dispersing low energy feature. An 
indication of this additional feature was also seen in the 
earlier low resolution data roughly 500 meV below the 
lowest energy feature whose origin was uncertain, but 
was attributed to a string resonance excitation^ Even 
though not as clear as in the Ca2Cu02Cl2 case, we also 
observed a similar high energy feature in Si^CuC^C^ 
(not shown). 

Due to the rapid dispersion this feature is most sys- 
tematically characterized by the Momentum Distribution 
Curves (MDCs) shown in figure 1(b). The MDC curves 
are well fit by a Lorentzian up to binding energies of 
2 eV at which point the main valence band dominates 
the spectral weight. The resulting peak positions deter- 
mined by fitting the MDCs are overlayed on figure 2 from 
which both the low and high energy dispersions discussed 
above are evident. One can also track the dispersion from 
the peak positions of the EDCs in figure 1(a) which are 
marked with circles and triangles. The dispersion ob- 
tained from the MDCs and EDCs between 0.3 and 0.7 
eV (low energy part) roughly agree. At high energies (~ 



1.3 eV) the MDCs also agree with the very broad and 
weak feature seen in the EDCs (triangles). The fact that 
the MDCs contain a single peak at all energies suggests 
that the suppressed spectral weight near the T point is 
not caused by the matrix element effect, but is due to 
a significant transfer of spectral weight from the low en- 
ergy feature to the rapidly dispersing high energy fea- 
ture. Also plotted in the figure is the dispersion from the 
LDA calculation shifted by ~ 0.7 eV (dashed line). Even 
though there is no independent justification for shifting 
the LDA dispersion by 0.7 eV, the agreement between 
the slopes of the two dispersions is still remarkable. This 
suggests that the rapidly dispersing high energy feature 
may be tracking the unrenormalizcd band. 

The low and high energy dispersive features form a 
kinklikc feature as can be seen in the image of figure 2. 
Ignoring the energy scale, this is very reminiscent of the 
data on metallic hole-doped cuprates, where a kink in 
the quasi-particle dispersion has been attributed to cou- 
pling to a phonon mode of ~ 70 meV^ In our data the 
kinklike feature occurs roughly 0.5 eV below the top of 
the valence band, which can be seen from the MDC de- 
rived dispersion in figure 2. However, in the present case 
the kink can not originate from coupling to a bosonic 
mode. First of all, the system under study is an insu- 
lator with a 1.5 eV charge transfer gap and coupling to 
a bosonic mode is expected to happen only in metal- 
lic systems in the standard theory. In addition, we also 
stress that, from the extremely large widths of the EDC's 
and MDCs, the features we are studying are certainly 
not quasi-particles^ In regards to the latter point, the 
MDC width of the high energy feature clearly decreases 
as the binding energy increases above 1 eV (see figure 
1(b)). This is counterintuitive as the phase space for 
decay channels of the photohole should increase as the 
binding energy increases, which should thus further raise 
a red flag when attempting to interpret these features as 
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FIG. 3: (a) EDCs along the r-(27r,0) direction taken with 
22.4 eV photons. The spectral weight at (27r,0) is weak, as is 
the case for F. (b) Image plot of ARPES data along the T- 
(0,7r) cut with 15.5 eV photons. The in-plane polarization of 
the light is along the (n,0) direction. Overlaid on the plot is 
the dispersion from LDA calculation shifted by 0.2 eV (dashed 
line) 

quasi-particle bands. 

The apparent transfer of spectral weight between the 
two observed dispersing features may give an answer to 
the puzzling observation of the vanishing spectral weight 
at the r point in all cuprates. Even though the dis- 
persions in cuprates, especially in insulating compounds, 
have been said to be renormalized to one with ~0.3 eV 
dispersion, the charge transfer peak mysteriously disap- 
pears near the T point. Traditionally, this lack of spec- 
tral weight at T has been attributed to matrix element 
effects. 5 The photoemission intensity is roughly expressed 
as / ~ |(i|A • p|/)| 2 where \i) and |/) are initial and fi- 
nal states, respectively. Both the initial \i) (primarily Cu 
3d x 2_ y 2) and final states |/} (« e lkz ) have even symmetry 
with respect to the (7r, 0) direction at the T point. Since 
A • p is a dipole operator, it provides an odd symmetry 
if the polarization is in the plane and along the (n, 0) 
direction. This has been the argument for the vanishing 
spectral weight at the T point. In real cases, however, 
such symmetry condition is never perfectly met, as the 
polarization almost always has an out-of-plane compo- 
nent (taken here as the z-direction)^ This additional 
^-component contributes to a non-vanishing matrix ele- 
ment, and thus the spectral weight should be finite at 
the T point. In fact, the manganites do possess signif- 
icant spectral weight at T due to a band with d 3z 2_ r 2 
orbital character for which the above arguments would 
be identical^ 



Perhaps even more importantly, we have failed to find 
any photoemission intensity at the equivalent fc-space lo- 
cation of (27T, 0) where the selection rule is no longer rel- 
evant. This can be seen in figure 3(a) where upon ap- 
proaching (2tt, 0) from (0, 0) spectral weight is rapidly 
pushed to higher binding energies, and the spectra at 
(2ir, 0) look identical to (0, 0). Furthermore, note that the 
matrix element argument dictates that the entire (0, 0) 
to (0, 7r) cut should have vanishing weight due to the fact 
that the in-plane component of polarization of the light 
is along the (it, 0) direction. Although it is indeed sup- 
pressed with respect to the (0, 0) to (n 7 0) cut, there still 
remains finite low energy weight above (0, 7r/4), as can 
be seen in figure 3(b). This not only vanishes upon ap- 
proaching the r point, but also has a high binding energy 
feature which disperses rapidly to higher binding energy. 
This is identical to the T to (7r, tt) cut, including the ob- 
servation that the rapidly dispersing high energy feature 
has a similar dispersion to that of LDA, as can be seen 
in panel 3(b). 

All of this suggests that in addition to the matrix el- 
ements which will suppress spectral intensity at the T 
point, there is a transfer of spectral weight to very high 
binding energies as the T point is approached. This is 
a necessary consequence of the ARPES sum rules which 
say that the total spectral weight of the single particle 
spectral function must be conserved (J A(k,uj)duj = 1). 
Note, this statement is true even for highly k dependent 
self energies. 

Theoretically, there are several potential explanations 
for the multiple features in the photoemission spectra. 
Cluster calculations identify the first hole ionization state 
as a singlet state and also calculate a large spectrum 
of excitations which one will observe in a photoemis- 
sion experiment. 17 While the majority of the states lie 
too high in energy (such as the Zhang-Rice triplet state 
at ~3.5 eV), there are states which are only ~1.5 eV 
above the lowest lying singlet state. Although we can 
not rule this out conclusively, these states are still some- 
what too large in energy, and there is also no a priori 
reason to expect these states to mimic the dispersion of 
the underlying band structure. Thus it is unlikely that 
they should correspond to our rapidly dispersing band- 
like feature. Alternatively, momentum resolved numer- 
ical calculations of the single band Hubbard model and 
the extended t-J model have large amounts of incoherent 
spectral weight (>50% at many fc-points) with significant 
structure in momentum space spread over large energies. 
Surprisingly, although approached differently from the 
numerical algorithm and approximations chosen, these 
calculated results bear similarities to what we observe 
in our data, whether it be from exact diagonalization, 18 
self-consistent Born approximation^ quantum Monte 
CarlojSS strong coupling perturbation theoryjSi or varia- 
tional cluster perturbation theory^ (a type of generalized 
dynamical mean field theory) . We also note that the high 
energy structure was sometimes referred to as "string" 
states in the earlier literature*^ The recent calculations 
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show as r is approached spectral weight is transferred 
from a low energy band that mimics that of the t-J 
model to spectral weight at higher energies with an en- 
ergy scale of the hopping energy t, similar to what we see 
experimentally. The two features have been interpreted 
in terms of the coherent motion of a hole in the antiferro- 
magnetic spin background and the incoherent motion of 
a hole hopping within the antiferromagnetic spin bag^S 
Thus, one might naturally expect the incoherent feature 
to track the unrenormalizcd band. 

The outstanding problem with most of the numerical 
calculations is that the lowest energy excitation still cor- 
responds to some sort of quasiparticle like pole, which is 
not consistent with ARPES experiments. The very large 
widths both in EDCs and MDCs are unphysical lifetimes 
for an associated quasiparticle. Also, the low-energy por- 
tion of the EDC lineshape can not be fit by any sort of 
power-law tail as would be expected from a quasiparticle 
pictured Finally, the doping dependence of the spec- 
tral weight argues that the quasi-particle states which 
are formed at the chemical potential in the hole-doped 
compounds receive spectral weight from the incoherent 
background which is the same feature we identify as the 
charge transfer band in the insulatoriSii2^ 

A fully consistent picture of the insulator must there- 
fore also be able to explain the large widths which are ob- 
served. An obvious candidate for the cause of such large 
widths is the strong electron-electron correlation, which 
suggests that a multiple state picture be naturally evoked 
in the photoemission spectra and that the true quasi- 
particle has vanishing weight^ In support of such a view, 
the phase string calculations based on the t-J model have 
a vanishing quasi-particle residue, and the spectra, com- 
posed entirely of incoherent weight, reasonably reproduce 
the experimental lineshapes and dispersion^ Alterna- 
tively, the multiple states that constitute the broad peaks 
observed in the experiments may be many shake-off peaks 
as in the case of the Franck-Condon effect^ The appeal- 
ing aspect of these types of pictures is that, as we have 
already demonstrated, our highly dispersing feature in 
Ca2Cu02Cl2 can most naturally be explained as being 
structure in the incoherent part of the spectral function. 
Indeed, in support of the latter view some calculations 



predict that the broad peaks still track the t-J disper- 
sion even when the polarons are of lattice origini^SiSS 
The large widths of the insulator can then be explained 
as resulting entirely from the structure in the incoherent 
spectral function which tracks both energy scales in the 
problem: a ' J'-like band at low energy and a T-like band 
at high energies. 

In conclusion, we use high resolution ARPES data of 
Ca2CuC>2Ci2 to elucidate a rapidly dispersing feature at 
high binding energies, which is very close to the unrenor- 
malized bare band. We find a transfer of spectral weight 
as a function of momentum from the low energy feature, 
which tracks the extended t-J model dispersion, to this 
high energy feature as T is approached. This interpreta- 
tion is consistent with several different calculations work- 
ing on the half-filled CuC>2 plane. This effect is in addi- 
tion to the matrix element effects which we believe are 
insufficient to explain the lack of spectral weight at T. 
The high-energy feature highlights the fact that the in- 
coherent spectral weight is spread over very large energies 
and has significant structure in momentum space. This 
potentially allows a previously overlooked aspect of the 
spectral function to be used to differentiate various mod- 
els which purport to understand the spectra of strongly 
correlated electron systems. 
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